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Abstract. Based on 3D hydrodynamical model atmospheres computed with the CO s BOLD 
code and 3D non-LTE (NLTE) line formation calculations, we study the effect of the convection- 
induced line asymmetry on the derived 6 Li abundance for a range in effective temperature, grav- 
ity, and metallicity covering the stars of the Asplund et al. (2006) sample. When the asymmetry 
effect is taken into account for this sample of stars, the resulting 6 Li/ 7 Li ratios are reduced by 
about 1.5% on average with respect to the isotopic ratios determined by Asplund et al. (2006). 
This purely theoretical correction diminishes the number of significant 6 Li detections from 9 
to 4 (2a criterion), or from 5 to 2 (3<r criterion). In view of this result the existence of a 6 Li 
plateau appears questionable. A careful reanalysis of individual objects by fitting the observed 
lithium 6707 A doublet both with 3D NLTE and ID LTE synthetic line profiles confirms that 
the inferred 6 Li abundance is systematically lower when using 3D NLTE instead of ID LTE line 
fitting. Nevertheless, halo stars with unquestionable 6 Li detection do exist even if analyzed in 
3D-NLTE, the most prominent example being HD 84937. 

Keywords, stars: abundances, stars: atmospheres, hydrodynamics, convection, radiative trans- 
fer, line: formation, line: profiles, stars: individual (G271-162, HD 74000, HD 84937) 



1. Introduction 

The spectroscopic signature of the presence of 6 Li in the atmospheres of metal-poor 
halo stars is a subtle extra depression in the red wing of the 7 Li doublet, which can only 
be detected in spectra of the highest quality. Based on high-resolution, high signal-to- 
noise VLT/UVES spectra of 24 bright metal-poor stars, Asplund et al. (2006) report the 
detection of 6 Li in nine of these objects. The average 6 Li/ Li isotopic ratio in the nine 
stars in which 6 Li has been detected is about 4% and is very similar in each of these stars, 
defining a 6 Li plateau at approximately logn( 6 Li) = 0.85 (on the scale logn(H) = 12). A 
convincing theoretical explanation of this new 6 Li plateau turned out to be problematic: 
the high abundances of 6 Li at the lowest metallicities cannot be explained by current 
models of galactic cosmic-ray production, even if the depletion of 6 Li during the pre- 
main-sequence phase is ignored (see reviews by e.g. Christlieb 2008, Cayrel et al. 2008, 
Prantzos 2010 [this volume] and references therein). 

A possible solution of the so-called 'second Lithium problem' was proposed by Cayrel 
et al. (2007), who point out that the intrinsic line asymmetry caused by convection in the 
photospheres of metal-poor turn-off stars is almost indistinguishable from the asymmetry 
produced by a weak 6 Li blend on a presumed symmetric 7 Li profile. As a consequence, the 
derived 6 Li abundance should be significantly reduced when the intrinsic line asymmetry 
in properly taken into account. Using 3D NLTE line formation calculations based on 3D 
hydrodynamical model atmospheres computed with the C0 5 BOLD code (Freytag et al. 
2002, Wedemeyer et al. 2004, see also |http : //www . astro . uu . se/~bf / co5bold_main . html ) , 
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we quantify the theoretical effect of the convection-induced line asymmetry on the re- 
sulting 6 Li abundance as a function of effective temperature, gravity, and metallicity, for 
a parameter range that covers the stars of the Asplund et al. (2006) sample. 

A careful reanalysis of individual objects is under way, in which we consider two al- 
ternative approaches for fixing the residual line broadening, Vbr, the combined effect of 
macroturbulence (ID only) and instrumental broadening, for given microturbulence (ID 
only) and rotational velocity: (i) treating Vbr as a free parameter when fitting the Li 
feature, (ii) deriving Vbr from additional unblended spectral lines with similar proper- 
ties as Li 1 6707. We show that method (ii) is potentially dangerous, because the inferred 
broadening parameter shows considerable line-to-line variations, and the resulting 6 Li 
abundance depends rather sensitively on the adopted value of Vbr. 



2. 3D hydro dynamical simulations and spectrum synthesis 

The hydrodynamical atmospheres used in the present study are part of the CIFIST 
3D model atmosphere grid (Ludwig et al. 2009). They have been obtained from real- 
istic numerical simulations with the C0 5 BOLD code which solves the time-dependent 
equations of compressible hydrodynamics in a constant gravity field together with the 
equations of non-local, frequency-dependent radiative transfer in a Cartesian box repre- 
sentative of a volume located at the stellar surface. The computational domain is periodic 
in x and y direction, has open top and bottom boundaries, and is resolved by typically 
140x140x150 grid cells. The vertical optical depth of the box varies from logTR OSS » —8 
(top) to logrR OSS « +7.5 (bottom), and the radiative transfer is solved in 6 or 12 opacity 
bins. Further information about the models used in the present study is compiled in Ta- 
bled Each of the models is represented by a number of snapshots, indicated in column 
(6), chosen from the full time sequence of the corresponding simulation. 

These representative snapshots are processed by the non-LTE code NLTE3D that 
solves the statistical equilibrium equations for a 17 level lithium atom with 34 line tran- 
sitions, fully taking into account the 3D thermal structure of the respective model atmo- 
sphere. The photo-ionizing radiation field is computed at 704 frequency points between 
A 925 and 32 407 A, using the opacity distribution functions oflCastelli & Kurucz (2004) 



to allow for metallicity-dependent line-blanketing, including the HI-H+ and HI-HI 
quasi-molecular absorption near A 1400 and 1600 A, respectively. Collisional ionization 
by neutral hydrogen via the charge transfer reaction H(ls) + Li(n£) <-> Li + (ls 2 ) + H~ is 



treated according to Barklem, et al. (2003) More details are given in Sbordone et al. (2009) 



Finally, 3D NLTE synthetic line profiles of the Li I A 6707 A doublet are computed with 
the line formation code Linfor3D (http: //www. aip.de/^mst/linf or3D_main.html), using 
the departure coefficients 6, = (NLTE) /n^ (LTE) provided by NLTE3D for each level i 
of the lithium model atom as a function of geometrical position within the 3D model 
atmospheres. As demonstrated in Fig.[TJ 3D NLTE effects are very important for the 
metal-poor dwarfs considered here: they strongly reduce the height range of line forma- 
tion such that the 3D NLTE equivalent width is smaller by roughly a factor 2 compared 
to 3D LTE. Ironically, the line strength predicted by standard ID mixing-length models 
in LTE are close to the results obtained from elaborate 3D NLTE calculations. We note 
that the half-width of the 3D NLTE line profile, FWHM(NLTE) = 8.5 km/s, is larger by 
about 10%: FWHM(LTE) = 7.7 and 7.5 km/s, respectively, before and after reducing the 
Li abundance such that 3D LTE and 3D NLTE equivalent widths agree. This is because 
3D LTE profile senses the higher photosphere where both thermal and hydrodynamical 
velocities are lower. However, the NLTE line profile is significantly less asymmetric than 
the LTE profile, even if the latter is broadened to the same half-width (Fig.Q] bottom 
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Table 1. List of models used in the present study. Columns (2)-(6) give effective temperature, 
surface gravity, metallicity, number of opacity bins used in the radiation hydrodynamics simu- 
lation, and number of snapshots selected for spectrum synthesis. The equivalent width of the 
synthetic 3D non- LTE 7 Li doublet at A 6707 A, assuming A(Li)=2.2 and no Li, is given in 
column (7). Columns (8) and (9) show g„(Li) and qJ(Li), the Li/ 7 Li isotopic ratio inferred 
from fitting this 3D non-LTE line profile with two different kinds of ID profiles, in each case 
assuming a rotational broadening of v sini = 0.0 / 2.0 km/s, respectively (see text for details). 



Model 


rp 1) 

-'off 

[K] 


logs 


[Fe/H] 


Nbin 


Nsnap 


W 2 > 
[mA] 


^=n( 6 Li)/n( 7 Li) 

(3D) NLTE [%] 


q* b = n( 6 Li)/n( 7 Li) 

ID LTE [%] 


d3t59g40mm30n02 
d3t59g45mm30n01 
d3t63g40mm30n01 
d3t63g40mm30n02 
d3t63g45mm30n01 
d3t65g40mm30n01 
d3t65g45mm30n01 


5846 
5924 
6269 
6242 
6272 
6408 
6556 


4.0 
4.5 
4.0 
4.0 
4.5 
4.0 
4.5 


-3.0 
-3.0 
-3.0 
-3.0 
-3.0 
-3.0 
-3.0 


6 
6 
6 
12 
6 
6 
6 


20 
19 
20 
20 
18 
20 
12 


44.9 
39.9 
23.9 
24.1 
24.3 
20.0 
16.4 


1.14 / 1.14 
0.75 / 0.75 
1.96 / 1.93 
1.81 / 1.80 
1.07 / 1.06 
1.75 / 1.70 
1.29 / 1.27 


0.88 / 0.88 
0.63 / 0.64 
1.86 / 1.83 
1.63 / 1.62 
1.02 / 1.00 
1.70 / 1.66 
1.25 / 1.22 


d3t59g35mm20n01 
d3t59g40mm20n02 
d3t59g45mm20n01 
d3t63g35mm20n01 
d3t63g40mm20n01 
d3t63g40mm20n02 
d3t63g45mm20n01 
d3t65g40mm20n01 
d3t65g45mm20n01 


5861 
5856 
5923 
6287 
6278 
6215 
6323 
6534 
6533 


3.5 
4.0 
4.5 
3.5 
4.0 
4.0 
4.5 
4.0 
4.5 


-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 


6 
6 
6 
6 
6 

12 
6 
6 
6 


20 
20 
18 
20 
16 
20 
19 
19 
19 


42.8 
45.2 
42.3 
22.1 
23.7 
25.1 
23.0 
16.2 
17.0 


2.48 / 2.44 
1.59 / 1.59 
1.12 / 1.13 
4.09 / 4.00 
1.95 / 1.94 
1.92 / 1.92 
1.18 / 1.18 
2.28 / 2.21 
1.31 / 1.29 


2.02 / 2.01 
0.96 / 0.98 
0.45 / 0.46 
4.04 / 3.97 
1.79 / 1.78 
1.66 / 1.67 
0.97 / 0.97 
2.22 / 2.16 
1.21 / 1.19 


d3t59g40mml0n02 
d3t59g45mml0n01 
d3t63g40mml0n01 
d3t63g40mml0n02 
d3t63g45mml0n01 
d3t65g40mml0n01 
d3t65g45mml0n01 


5850 
5923 
6261 
6236 
6238 
6503 
6456 


4.0 
4.5 
4.0 
4.0 
4.5 
4.0 
4.5 


-1.0 
-1.0 
-1.0 
-1.0 
-1.0 
-1.0 
-1.0 


6 
6 
6 

12 
6 
6 
6 


20 
08 
20 
20 
20 
20 
19 


41.5 
38.2 
22.0 
21.8 
23.4 
15.5 
17.1 


1.60 / 1.62 

1.14 / 1.15 
2.37 / 2.38 

2.15 / 2.16 
1.36 / 1.37 
3.10 / 3.02 
1.51 / 1.51 


1.45 / 1.47 
0.83 / 0.85 
2.33 / 2.33 
2.05 / 2.06 
1.23 / 1.24 
3.14 / 3.06 
1.44 / 1.43 



Notes: 1 - ) averaged over selected snapshots; 2 ' averaged over selected 3D non-LTE spectra 
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3. 6 Li bias due to convective line asymmetry 

As outlined above, the 6 Li abundance is systematically overestimated if one ignores the 
intrinsic asymmetry of the 7 Li line components. To quantify this bias theoretically, we rely 
on synthetic spectra. The idea is as follows: we represent the observation by the synthetic 
3D NLTE line profile of the Li line blend, computed with zero 6 Li content. Except for 
an optional rotational broadening, the only source of non-thermal line broadening is 
the 3D hydrodynamical velocity field, which also gives rise to a convective blue-shift 
and an intrinsic line asymmetry. Now this 3D 7 Li line blend is fitted by 'classical' ID 
synthetic line profiles composed of intrinsically symmetric components of 6 Li and 7 Li. 
Four parameters are varied independently to find the best fit (minimum x 2 )- m addition 
to the total 6 Li+ 7 Li abundance, A(Li), and the 6 Li/ 7 Li isotopic ratio, <?(Li), which control 
line strength and line asymmetry, respectively, we also allow for a residual line broadening 
described by a Gaussian kernel with half-width Vbr, and a global line shift, Av. Note that 
the four fitting parameters are non-degenerate, since each one has a distinctly different 
effect on the line profile. The rotational line broadening is fixed to the value used in the 
3D spectrum synthesis (we tried v sini = and 2 km/s). The value g* (Li) of the best fit 
is then identified with the correction that has to be subtracted from the 6 Li/ 7 Li isotopic 
ratio derived from the ID analysis to correct for the bias introduced by neglecting the 
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Figure 1. Comparison of 3D LTE (dashed), 3D non-LTE (thick solid), and ID LTE (thin 
solid) equivalent width contribution functions (top) and line bisectors (bottom) of a sing le 7 Li 
component, computed for a typical metal-poor turn-off halo star (T e g = 6215 K, log g — 4.0, 
[Fe/H] = —2). Non-LTE effects strongly reduce the height range of formation and equivalent 
width of the line (W: 35.5 -)■ 15.6 mA), while differences between ID LTE and 3D NLTE are 
much smaller. The line asymmetry is smaller in non-LTE: the velocity span of the bisector 
is <5-u(NLTE) w 190 m/s, while <5u(LTE) « 310 and 260 m/s, respectively, before and after 
adjusting the equivalent width (by reducing the Li abundance) and the half-width (by Gaussian 
broadening) of the LTE line to match the values of the NLTE profile (compare the two dashed 
bisectors). 
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Figure 2. 6 Li/ 7 Li iso topic ratio, and ±lcr error bars, as a function of metallicity as derived by 
Asplund et al. (2006) before (top) and after (bottom) subtraction of <j£(Li) (Table[T] col. (9)) 



to correct for the bias due to the intrinsic line asymmetry. Open circles denote non-detections, 
filled circles and filled circles with a black dot denote 6 Li detections above the 2a and 3a level, 
respectively. 



intrinsic line asymmetry: q( 3D ^(Li) = q( 1D \~Li) - q*(Li). The procedure properly accounts 
for radiative transfer in the lines, including saturation effects. 

Two different sets of ID profiles were used for this purpose: (a) NLTE line profiles 
based on the (3D) model, constructed by averaging the 3D model on surfaces of constant 
optical depth, and (b) LTE line profiles computed from a so-called LHD model, a ID 
mixing-length model atmosphere that has the same stellar parameters and uses the same 
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microphysics and radiative transfer scheme as the corresponding 3D model. For both 
kind of ID models, the microturbulence was fixed at £ m ; c = 1.5 km/s. The mixing length 
parameter adopted for the LHD models is «mlt = 0.5. 

We have determined q* (Li) according to the method outlined above for our grid of 3D 
model atmospheres. The results are given in Table[TJ both for fitting with the (3D) NLTE 
profiles (q*(Li), col. (8)) and with the ID LTE profiles of the LHD models (g£(Li), col. (9)). 
At given metallicity, the corrections are largest for low gravity and high effective temper- 
ature. They increase towards higher metallicity. We also note that q* (Li) is essentially 
insensitive to the choice of the rotational broadening. 

The analysis of Asplund et al. (2006) utilizes ID LTE profiles computed from MARCS 
model atmospheres. Hence, the correction g^(Li), computed with the ID LTE profiles 
of the ID LHD models, should be applied to their 6 Li/ 7 Li isotopic ratios. The resulting 
downward corrections are typically in the range 1% < gJ(Li) < 2% for the stars of their 
sample (cf. Steffen et al. 2010, Fig. 2). After subtracting the individual <j£(Li) for each of 
these stars, according to their T e g, logg, and [Fe/H], the mean 6 Li/ 7 Li isotopic ratio of 
the sample is reduced from 0.0212 to 0.0056, as illustrated in Fig. [2] Given the 6 Li/ 7 Li 
isotopic ratios and their la error bars as determined by Asplund et al. (2006)| the num- 
ber of stars with a 6 Li detection above the 2a and 3a level is 9 and 5, respectively, out of 
24. After correction, the number of 2a and 3a detections is reduced to 4 and 2, respec- 
tively. One of the stars meeting the 3a criterion after correction, HD 106038, survives 
only because of its particularly small error bar of ±0.006. The other one is G020-024, 
which shows the clearest evidence for the presence of 6 Li (q(Li) = 0.052 ± 0.017), while 
HD 102200 remains the only clear 2a detection (g(Li) = 0.033 ± 0.013). The spectra of 
these stars should be reanalyzed with 3D NLTE line profiles. 



4. Analysis of observed spectra: 3D NLTE versus ID LTE line fitting 

As a further exercise, we have fitted the observed Li I A 6707 A spectra of three halo 
turn-off stars (see Tabled, both with 3D NLTE and ID LTE synthetic line profiles. Since 
the parameters of these three stars are very similar, the synthetic spectra were computed 
in all cases from hydrodynamical model d3t6300g40mm20n01 (see Table[T]) and the cor- 
responding ID LHD model. The standard approach is to vary the four parameters A(Li), 
g(Li), Vbr, and Av to find the optimum line profile fit, as described above (method I). 
The results are presented in TableO As expected, the 3D analysis yields lower 6 Li/ 7 Li 
isotopic ratios by about 1.7%. The numbers differ slightly from those given by Steffen 
et al. 2010 due to an upgrade of the line fitting procedure, but the conclusions remain 
unchanged: HD 74000 and G271 — 162 are considered non-detections, while HD 84937 re- 
mains a clear 6 Li detection with q( 3D ~ NLTE) (Li) » 5%. 

For completeness, Table[3] also shows the results of fitting the observed lines with 3D 
LTE synthetic profiles (not recommended). The 6 Li abundances obtained in this way lie 
between the 3D NLTE and the ID LTE results. It is not obvious why fitting with the 
more asymmetric 3D LTE profiles gives higher 6 Li than fitting with the slightly more 
symmetric 3D NLTE profiles, but this is a robust result. 

For HD 74000, we also tried an alternative approach (method II), where Vbr is deter- 



Table 2. Observed Li I A 6707 A spectra of three metal-poor turn-off stars 



Star 


T cff [K] 


log9 


[Fe/H] | R=A/AA 


S/N 


Instrument 


Reference 


HD 74000 
G271-162 
HD 84937 


6203 
6330 
6300 


4.03 
4.00 
4.00 


-2.05 
-2.25 
-2.30 


120 000 
110 000 
100 000 


600 
600 
630 


ES03.6 / HARPS 

VLT / UVES 
CFHT / GECKO 


Cayrel et al. 2007 
Nissen et al. 2000 
Cayrel et al. 1999 
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Table 3. Fitting the observed spectra of Table[2]with 3D and ID synthetic line profiles. 
Columns (4)-(7) show the results for v sin i = 0.0 / 2.0 km/s. 



Star 


synthetic 


fitting 


A(Li) 1} 


n( 6 Li)/n 


( 7 Li) 


V 2) 

BR. 


Av 




spectrum 


method 




[%] 




[km/s] 


[km/s] 


HD 74000 


3D NLTE 


I 


2.25 / 2.25 


-1.1 / 


-1.1 


3.1 / 2.1 


0.64 / 0.64 




(3D LTE) 


I 


1.85 / 1.85 


-0.4 / 


-0.5 


4.7 / 4.1 


0.66 / 0.67 




ID LTE 


I 


2.23 / 2.23 


0.6 / 


0.6 


5.9 / 5.4 


0.43 / 0.43 




3D NLTE 


II 


2.25 / 2.25 


-1.1 / 


-1.3 


3.1 / 2.2 


0.64 / 0.65 




(3D LTE) 


II 


1.84 / 1.84 


4.6 / 


4.2 


3.1 / 2.2 


0.46 / 0.47 




ID LTE 


II 


2.23 / 2.23 


1.7 / 


1.7 


5.6 / 5.1 


0.39 / 0.39 


G271-162 


3D NLTE 


I 


2.30 / 2.30 


0.7 / 


0.6 


3.7 / 2.9 


0.04 / 0.04 




(3D LTE) 


I 


1.89 / 1.89 


1.3 / 


1.1 


5.2 / 4.6 


0.06 / 0.07 




ID LTE 


I 


2.27 / 2.27 


2.3 / 


2.4 


6.2 / 5.8 


-0.17 / -0.17 


HD 84937 


3D NLTE 


I 


2.20 / 2.20 


5.2 / 


5.2 


3.3 / 2.4 


0.02 / 0.02 




(3D LTE) 


I 


1.80 / 1.80 


5.8 / 


5.7 


4.8 / 4.3 


0.05 / 0.05 




ID LTE 


I 


2.18 / 2.18 


6.9 / 


6.9 


6.0 / 5.6 


-0.19 / -0.19 



Notes: 1} log [n( 6 Li) + n( 7 Li)] - logn(H) + 12: 2) Gaussian kernel, bold: fixed from Fel lines 

mined from other spectral lines, and is then fixed during the fitting of the Li I A 6707 A 
feature. For this purpose, we rely on a set of 6 clean Fe I lines with similar excitation 
potential between 2.4 and 2.6 eV: the 5 lines of Cayrel et al. (2007), Table 1, plus Fel 
A 6230.7 A, Ei = 2.559 eV. Their equivalent widths range from W ~ 13 to 30 mA, em- 
bracing the strength of the Li doublet. Fitting the 6 Fel lines with 3D LTE synthetic 
line profiles, we obtain VeR(Fe) = 3.05 (2.15) ±0.16 (0.23) km/s for v sini = (2) km/s, 
in close agreement with Vbr(LI) inferred from the fitting of the Li I profile. This result 
may be taken as an indication that, in contrast to Li, the selected Fel lines are not 
severely affected by departures from LTE. We note, however, a trend of VeR(Fe) in- 
creasing slightly with W, which remains to be understood. If instead the Fel lines are 
fitted with ID LTE synthetic profiles, we find VeR(Fe) = 5.61 (5.13) ±0.07 (0.09) km/s, 
systematically lower than Vbr(LI) by 0.3 km/s. Fixing VBR=VBR,(Fe), the best fit of 
the Lil doublet implies q 1D LTE (Li)=1.7%, which is 1.1% higher than with method I. 
Finally, Table|3] demonstrates that applying method II with 3D LTE fitting of Li I leads 
to a severe overestimation of the 6 Li/ 7 Li isotopic ratio: q 3D LTE (Li) > 4% compared to 
q 3B NLTE (Li)=-l.l%! 



5. Conclusions 

The 6 Li/ 7 Li isotopic ratio derived from fitting of the Lil doublet with 3D NLTE syn- 
thetic line profiles is shown to be about 1% to 2% lower than what is obtained with ID 
LTE profiles. Based on this result, we conclude that only 2 out of the 24 stars of the 



Asplund et al. (2006) sample would remain significant 6 Li detections when subjected to 
a 3D non-LTE analysis, suggesting that the presence of 6 Li in the atmospheres of galactic 
halo stars is rather the exception than the rule, and hence does not necessarily constitute 
a cosmological 6 Li problem. If we adopt the approach by Asplund et al. (2006)| relying 
on additional spectral lines to fix the residual line broadening, the difference between 3D 
NLTE and ID LTE results increases even more, as far as we can judge from our case 
study HD 74000. Until the 3D NLTE effects are fully understood for all involved lines, 
we consider this method as potentially dangerous. 
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